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A procedure to derive the Ginzburg-Landau (GL) theory from the multiband BCS Hamihonian 
is developed in a general case with an arbitrary number of bands and arbitrary interaction matrix. 
It combines the standard Gor'kov truncation and a subsequent reconstruction in order to match 
accuracies of the obtained terms. This reconstruction recovers the phenomenological GL theory as 
obtained from the Landau model of phase transitions but offers explicit microscopic expressions for 
the relevant parameters. Detailed calculations are presented for a three-band system treated as a 
prototype multiband superconductor. It is demonstrated that the symmetry in the coupling matrix 
may lead to the chiral ground state with the phase frustration, typical for systems with broken 
time-reversal symmetry. 
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I. INTRODUCTION 

Studies of multiband (or multigap) superconducting 
systems, where more then one carrier band contributes 
to the formation of the condensate, have now more then 
a half-century historyii"— In the last decade clear ex- 
perimental evidences of multigap condensates were ob- 
served in a rich variety of materials such as magne- 
sium diboridc)^ oxypnictidcsj^ iron arsenides^ and iron 
pnictidcsi^ The string of discoveries continues today so 
that the number of multiband superconductors increases 
almost yearly. 

There are different physical mechanisms responsible for 
the formation of multiple carrier bands. For example, in 
bulk specimens the multi-gap structure can be related 
to the appearance of separate pockets in the Fermi sur- 
face centered around some points of the Brillouin zone. 
However, it was recently shown that multiband super- 
conductivity can also arise in nanoscalc specimens (e.g. 
in nanofilms) made of ordinary single-band supercon- 
ducting materials, where the geometrical size quantiza- 
tion creates distinct carrier subbandsi^ These and similar 
findings broadened the interest in the physics of multi- 
band coherent phenomena, and that interest has given a 
strong impetus to theoretical investigations. One of the 
focal points of such investigations is revisiting long estab- 
lished and widely used theoretical models and methods in 
superconductivity in the context of multiband supercon- 
ductors. One of such methods is the Ginzburg-Landau 
(GL) theory^ which is commonly regarded as one of the 
most general and yet simple approaches for conventional 
single-band superconductors. Surprisingly, the general- 
ization of the GL theory to the multiband case is still a 
highly debated issue. 

On the microscopic level a multiband superconductor 
is modeled by the multiband generalization of the BCS 
theoryi^ The corresponding GL equations are derived 
using a straightforward application of the original single- 



band recipe by Gor'kov jii see, e.g., Refs. Il2lll3l . In this 
approach superconducting gap functions in each carrier 
band, hereafter referred to as band gaps, are regarded as 
the order parameters. Similarly to the single-band case, 
the anomalous Green's function of each band is expanded 
in powers of the corresponding band gap and its spatial 
gradients fii and then the expansion is truncated to keep 
the same terms as in the single-band GL theory. This 
procedure yields a system of nonlinear GL-likc equations, 
one for each band gap, coupled via the linear Josephson- 
like terms, and the corresponding multi-component func- 
tional. This is often referred to as the multi- component 
GL model and is widely used in the analysis of multiband 
superconductors ji^Ti^^ 

Although this formulation of the GL theory appears 
intuitively justified, partially by a familiar structure of 
the obtained equations, it possesses several fundamental 
inconsistencies. First, it has to be reconciled with the 
phenomenological Landau theory of phase transitions, 
according to which the order parameter must be asso- 
ciated with a particular irreducible representation of the 
relevant symmetry group. Following this prescription, 
Volovik and Gor'kov developed a classification of the ex- 
otic superconducting phases within the GL theory^^ (a 
systematic classification of the GL theories based on the 
symmetry analysis can be found in Ref. |24[ ) . It is impor- 
tant that the number of independent order parameters 
in the GL theory, given by the dimensionality of the irre- 
ducible representation, is typically lower then the number 
of bands, which is certainly different from the multicom- 
ponent model mentioned above. 

Second, the analysis of the multi-component GL model 
presented by Geilikman, Zaitscv and Kresin^ and more 
recently by Kogan and Schmalian;^ revealed another in- 
consistency: the accuracy of a solution to the formalism 
exceeds the accuracy of its derivation. This discrepancy is 
intrinsic in the multiband generalization of the Gor'kov 
procedure and can only be eliminated by invoking an 
additional truncating reconstruction, which removes the 



artificial higher-order contributions <^i^i^ — 

Without additional symmetries such reconstruction 
yields a strict proportionality of all band gaps, i.e., the 
GL theory has a single order parameter. This conclu- 
sion agrees with the phenomenological classification that 
predicts a single-component GL theory in this case. (De- 
viations from this result appear only in higher-order cor- 
rections to the ordinary GL theory»^"— ) 

This analysis did not consider the case of a degenerate 
solution for T^ which appears due to an additional sym- 
metry of the system. Furthermore, the calculations in 
Ref . y for an arbitrary number of bands employed a rather 
restrictive ansatz for the band gaps, while Refs. l26ll27l uti- 
lized the separability specific to the two-band case. The 
microscopic derivation of the multiband GL theory has 
not been yet achieved in the general case. Notice, that 
a mechanical merge of the symmetry analysis with the 
Gor'kov truncating procedure, in which the outcome of 
the Gor'kov procedure is simply rewritten in terms of the 
basis states of the relevant symmetry group representa- 
tions, does not solve the problem. It yields a mixture 
of different irreducible representations, which should not 
happen in the standard GL formalismi^ 

In this work we derive the reconstructed (true) GL the- 
ory from the microscopic Hamiltonian for a multiband su- 
perconductor in a general case with an arbitrary number 
of bands as well as with an arbitrary symmetry (reflected 
in the degeneracy of the solution for Tc). The origin of 
the symmetry is not important here. We note that it can 
appear not only due to the lattice structure of the mate- 
rial, as discussed in Ref. [23ll24l but also due to other rea- 
sons, e.g., the geometrical shape of the sample like in su- 
perconducting single-crystalline nanofilms.— A detailed 
analysis of the obtained equations is then performed for 
the three-band system treated as a prototype of a multi- 
band superconductor. In particular, we consider a simple 
three-band model of pnictides with dominant interband 
couplings which allows for the two-fold degeneracy of the 
solution for Tc- We demonstrate that in full agreement 
with the phenomenological GL theory this system has 
two order parameters, related to the two-dimensional ir- 
reducible representation of the relevant symmetry group. 
However, unlike the phenomenological analysis based on 
the symmetry consideration, the derivation from the mi- 
croscopic theory offers the explicit expressions for the co- 
efhcients of the GL theory. These expressions are highly 
nontrivial in the case of multiband superconductors be- 
cause they contain important information about contri- 
butions of different bands that cannot be obtained from 
the symmetry arguments j^ The corresponding ground 
state of the system is found to be a chiral state with a 
nontrivial phase difference between the band gaps. Such 
states in multigap superconductors have attracted much 
interestiiSi"— i^ as they could lead to unconventional phe- 
nomena such as the formation of antiferromagnetic do- 
mains or noninteger vortices, see, e.g., Ref. [2J. Notice 
that the present work does not go beyond the standard 
GL domain (i.e., band gaps are proportional to t^''^ , with 



T = 1 — T/Tc the proximity to the critical temperature) . 
An extended version of the multiband GL theory with 
the proper higher-order contributions to the band gaps 
will be published elsewhere. 

The paper is organized as follows. In Sec. |ll] the GL 
theory for multiband superconductors is derived starting 
from the standard multiband BCS model. The deriva- 
tion is performed in three steps: (i) the truncated multi- 
band gap equation is obtained in matrix form by fol- 
lowing the Gor'kov procedure adapted for the case of 
multiple bands, (ii) the truncating reconstruction is then 
applied by invoking the r-expansion, and (iii) an explicit 
form of the resulting GL equations is obtained by keep- 
ing the terms of order r^/^ in the band gaps. In Sec. IIIII 
we recast the final formalism in a more explicit form, for 
both the nondcgencratc and degenerate cases. In Sec. IIVI 
we consider a three-band model, for which expressions 
for the coefficients of the GL equations can be calculated 
analytically for an arbitrary interaction matrix. Then, 
we investigate the case of a degenerate solution for Tc 
for a simple variant of the model with strong interband 
couplings and demonstrate analytically that the degen- 
eracy in this model leads to the chiral ground state. Our 
summary and conclusions can be found in Sec. |Vl 



II. DERIVATION OF THE GL THEORY 



A. Truncated gap equation 

Following Gor'kov^ the GL theory is usually derived 
from the gap self-consistency equation, using an expan- 
sion of the anomalous Green's function in powers of the 
order parameter and its spatial derivatives. We out- 
line this derivation for multiple bands, starting from the 
multiband BCS Hamiltoniaiii'^ with the s-wave singlet 
pairing, which reads as 

Hbcs ^Hc + ^ jd\ [^ i^liv) r,(r)V;..(r) 



+ V'lt(r)'Ali(r)A,(r)-|-h.c. 



(1) 



where i enumerates different bands, il)ia{v) and ^'J^(r) are 
the carrier field operators, Ai(r) are the band gap func- 
tions, or simply band gaps. He is the c-number term (see, 
e.g., Ref. [l3l) . and T'i(r) stands for the single-electron en- 
ergy. Equation ([T]) is accompanied by the self-consistency 
gap equation 



A = ffi?, 



(2) 



where we introduce 5, the matrix of the coupling 
constants gij, and use the vector notations A"^ = 
(Ai,A2,...) and R^ = {Ri,R2,...), with i?, == 
{tljif{r)tpii{r)) being the anomalous Green's function of 
the z-th bandit 

Using Eq. ([T]), we expand the anomalous Green's func- 
tions in the vicinity of Tc into a series in powers of the 



band gaps and their spatial gradients. As the Hamilto- 
nian ([1]) is diagonal over the band index, the series is 
obtained independently for each band yielding the same 
expressions as in the single-band case. Referring inter- 
ested readers to the original calculations fii here we quote 
the well-known final expansion for Ri, where only the 
leading nonlinear and gradient terms arc retained, 



R,[A,]c^N,{0)AA,+n,[A,], 



where Ni{0) is the band DOS, A = ln(2^^), and 
n,[A,] = -a,A, - b,\A,\^A, + JC^V^A,. 



(3) 



(4) 



In Eqs. (|31) and ^ Uc is the cut-off frequency, F — 0.577 
is the Euler constant and the coefficients are calculated 



a, = -iV,;(0)T, h 






(5) 



with C(- • ■ ) the Riemann zeta function and Vi the band 
Fermi velocity. Although here only results for the clean 
limit are quoted, we note that the structure of the equa- 
tions will be the same for dirty systems, as is usually the 
case in the standard GL formalism.— 

We note that the magnetic field is not included in 
Eq. ((4]). The generalization to the nonzero-field case is 
trivial and will be done on a later stage by using the stan- 
dard prescription of inserting the gauge invariant gradi- 
ent D = V — jffA, where A is the vector potential. One 
should remember, however, that this recipe is valid ex- 
clusively for the standard GL domain when only terms of 
order r^ " are kept in the gap functions. A more involved 
and complex procedure is needed when higher-order cor- 
rections to the band gaps are incorporated.— 

Substituting Eq. ([3]) into Eq. ([2]), we obtain the follow- 
ing system of coupled equations: 

(A + a,) A, + 6^A,|A,p - /C^V^A, + ^ 7„Aj- = 0, 

(6) 

where 7y is the element of the inverted interaction matrix 
g~^ and the constants At are defined as 



A., = -/^r - N^{0)A. 



(7) 



The truncated equations given by Eg. ([6]) are com- 
monly referred to as a generalization of the GL theory to 
the multiband case, or the multi- component GL theory. 
This interpretation is suggestive, especially given that in 
the limit of zero interband couplings Eq. (j6]) yields N 
uncoupled GL equations for A,;'s that are the true Lan- 
dau order parameters for the uncoupled bands. Thus, the 
coupling is commonly assumed to be a weak perturbation 
that does not significantly alter this physical picture. 

However, already this trivial limit highlights shortcom- 
ings of the interpretation of Eq. ([H]) as a consistent multi- 
band GL formalism. In the absence of coupling each band 



has its own critical temperature Td while Tc of the en- 
tire system is the largest of those. In the vicinity of Tc, 
which is the usual validity domain of the GL theory, only 
the band with Td = Tc develops a superconductive state 
and, therefore, the system is in fact described by a single- 
band GL theory with the single order parameter (here we 
assume that Tci's are well separated). One can also imag- 
ine a degenerate situation when M < N gaps have the 
same largest critical temperature Tc. Here, in the zero- 
coupling limit, the system is described by the theory with 
M < N order parameters corresponding to M uncoupled 
components. Thus, in the zero interband-coupling limit, 
the GL theory always has fewer active order parameters 
then the number of the available bands. This conclusion 
is of course trivial in the nonintcracting case. However, 
in what follows we demonstrate that it holds also in the 
general case of a nonzero coupling within the accuracy of 
the GL approach. 



B. Reconstructed GL theory 

Deriving the GL theory for the general case of nonzero 
interband interactions starts by noting that as discussed 
in the introduction, Eq. ([6]) is inconsistent because the 
accuracy of its solution exceeds the accuracy of its 
derivation. One can see this (details can be found in 
Refs. I3ll4l l26l - l28r ) by taking into account that the coeffi- 
cients Ai + Oi and 7ij are not zero in the limit r — > 
(T -^ Tc). This implies that a solution to Eq. ([5]), 
when being expanded in r, comprises terms of arbitrar- 
ily high orders, i.e., all A^'s are given by infinite series 
in powers t""*"^/^, with integer n. At the same time the 
Gor'kov truncation neglects terms that contribute to or- 
ders higher than t^^^ in the band gap functions. The 
only situation when this inconsistency does not happen 
is the single-band GL theory where the coefficient of the 
linear term in the GL equation is proportional to t and, 
as a result, the solution comprises a single contribution 
of order t^/^. 

In order to reconcile the accuracy of a solution for A^ 
with the accuracy of the derivation of Eq. ^, we use 
the reconstruction procedure that abandons incomplete 
higher-order contributions from the band gaps. This pro- 
cedure is nothing more than a systematic perturbation 
expansion in r, which gives the GL theory and its cor- 
rections in a systematic wayi^ Following this procedure, 
the solution to Eq. ([6]) is sought in the form of a series in 
odd powers of r^/^ as 



A, = Af)+Af)+0(r5/2), 



(8) 



.(0) 



where A^^'^'' ex t^^"^ and A^ ' « r*^/^. This series is in- 
serted into Eq. ([5]) and then the terms of the same or- 
der are collected. A simple power-counting shows that 
making a solution to Eq. ([6]) consistent with the Gor'kov 
truncation, one should keep the two lowest orders in the 
resulting r-expansion of Eq. ^ and the leading order 
term ex r^/^ jj^ gq (p^ 



Notice that one must also take into account that spa- 
tial derivatives of the difference SAi = Ai — A^- do not 
contribute to these lowest orders. In an earlier consider- 
ation^ it was assumed that SAi is independent of the co- 
ordinates, see Eq. (14) in Ref. |j. Subsequent worka^l"— 
have demonstrated that such a restrictive ansatz is not 
needed. The GL theory introduces the coherence length, 
5 oc r~^/^ so that all spatial derivatives of all con- 
tributing terms in the band gaps scale as ex r^". In 
other words, each gradient operator V introduces a fac- 
tor ex T^'^. Counting powers of t in the expansion con- 
firms that the higher-order gradients of A^ as well as 
the lowest gradients of SAi do not contribute into the 
two lowest orders of the r-expansion of Eq. (jH]). 

Substituting Eq. ^ into Eq. ([6|) and collecting the 
terms of order r^' ^ we obtain the first equation in the 
reconstructed theory 



LA(") = 0, 



(9) 



where elements of matrix L are written as 



Lij = S^jAt + (1 - Sij)^ij (10) 



with Sij being the Kronecker symbol. The condition of 
solvability of Eq. dH), 



dctL = 0, 



(11) 



is the equation for the critical temperature Tc that gen- 
erally has N solutions. Clearly, one has to choose the 
solution with the maximal Tc as it yields the minimal 
value of the free energy. Equation (jH]) is commonly re- 
ferred to as the linearized gap equation as it can also be 
obtained by simply neglecting all the nonlinear contribu- 
tions in Eq. ([6]). 

When N > 2, one may encounter a situation with M < 
N degenerate solutions to Eq. pT|) that correspond to the 
same maximal value of Tc . In this case the matrix L has 
M eigenvectors ^q, with a = l,...,Af, corresponding to 
the zero eigenvalue of L at T = Tc. Without loss of 
generality these eigenvectors can be chosen orthogonal, 
and their normalization is not important. 

A general solution to Eq. © is then represented as a 
sum of M terms {M = 1 is for the nondegenerate case), 
one for each eigenvector, as 



AW(r)=^^„(r)e; 



(12) 



Here M functions V^q are specified by the equation that 
is obtained from Eq. (jS]) by matching terms of order r^' ^ 
as 

LA(i) =f1[AW], (13) 

where the components ^i[Al ] of J7[A^*'^] are given by 
Eq. (gl) with Ai replaced by A^^ A closed set of M 



equations for V'Q(r) is derived by projecting Eq. ((T3)) to 
the eigenvectors ^q,, which yields M equations given by 



j2u^A^?'>]^o, 



(14) 



where ^ai is the i-th component of ^^ . 



III. EXPLICIT FORM OF THE GL EQUATIONS 



A. Nondegenerate case 



Here we recast Eq. (|14p in a more explicit and famil- 
iar form. In the nondegenerate case a single function 
■(/)(r) = ipi{r) controls the same spatial profile of all band 
condensates. Rewriting Eq. (fT4|) for ip{r) one obtains 



ai' + b\i}\'^i;-IC'D^i; = 0, 



(15) 



where we include a nonzero magnetic field by replacing 
V — >■ D. The coefficients a, b and K. in Eq. p5|) are given 
by 



^a,C',& = E^»^-^ = E^^^^ 



(16) 



where ^i is the band component of ^ = ^i. The corre- 
sponding free-energy functional reads as 



J" 



d-^r 



a\^\' 



M' + IC\T)i;\' 



5! 



(17) 



Using this functional one derives the accompanying 
Maxwell equation for the gauge field in the form 

^rotB = i|^/C(^D>*-V*DV'). (18) 

As seen, Eqs. ([15]), ^7^ and ([TH]) have the form of the 
ordinary single-band GL theory. In fact, however, this 
is an effectively single-band GL theory as the coefficients 
a, b and JC comprise contributions of all bands. It is 
also important to remember that ip itself cannot be in- 
terpreted as an excitation gap: it is related to the band 
gap functions via Eq. P^ . 

The single-band representation of the reconstructed 
GL theory allows one to define the characteristic lengths 
of a multiband superconductor in a unique way. In par- 
ticular, the coherence length ^, the magnetic penetration 
depth A and their ratio k are given by the standard GL 
expressions as 




A = 



327r/C|ar 



K = $r 



327r3/C2 



(19) 



where $o is the fiux quantum. However, the multiband 
origin of Eqs. P^ - ((T5)) is still rcficcted in some properties 
of the system. For example, following Eqs. p^ and (|19p . 



one concludes that b, K, can be roughly estimated as lin- 
early proportional to A'^. Taking into account the relation 
K ex yJhjK? one arrives at the trend k ^ 1/^/N, which 
means that a multiband superconductor should approach 
the type-I character when the number of bands is large 
enough. 



B. Degenerate case 

When the maximal solution to Eq. pT|) is degenerate, 
i.e., M > 1, Eq. P^ yields a set of coupled nonlinear 
equations, an explicit form of which is obtained as 

y^ {aaP - ICap'D'^)lpp + ^ bal3^s 1pp1p*^1p5 = 0, (20) 

fs MS 

where the coefficients are defined as 

i i 

i 

The corresponding free-energy functional is now ob- 
tained in the form 



J": 



aP 

1 B^ 

-^ Y ba^S C V'/3 '07 ^5 + -^ 



oi^'yS 



(22) 



By calculating the functional derivative with respect to 
the vector potential, we obtain from Eq. ([22| the accom- 
panying Maxwell equation as 

^rotB = z^^/C„^(^„D*V;-V^DV'a). (23) 

ap 

The number of components in the reconstructed GL 
theory is 1 < M < TV, unlike in the original system of 
equations given by Eq. ^. Another important differ- 
ence is that all coefficients of the linear terms in the re- 
constructed GL theory are now proportional to r, which 
dictates that -00 oc r^/^. This eliminates the problem of 
the mismatch between the accuracy of the solutions and 
equations, which was the reason to seek the reconstruc- 
tion. 

Notice that the reconstructed GL formalism, obtained 
here by the r-expansion, recovers the standard Landau 
theory of phase transitions. In particular, the degener- 
acy of the linearized gap equation is related to an extra 
symmetry between bands, hidden in the relevant cou- 
pling matrix. The degree of degeneracy M is defined by 
the dimensionality of the corresponding irreducible rep- 
resentation with the M basis vectors ^q's. Equation ([^^ 
is interpreted as the Landau free-energy functional with 
tpaS being Landau order parameters. The reconstruction 



can thus be regarded as the procedure of finding the true 
Landau order parameter of the system, in the form of 
linear combinations of the band gaps, see Eq. p^ . How- 
ever, Eqs. (|17p and (|22p are derived by matching all rel- 
evant terms in the r-expansion, rather than through the 
phcnomenological approach based on the group-theory 
analysis)^ 

In agreement with the Landau recipe, the recon- 
structed GL theory is based on a single irreducible repre- 
sentation. However, if one continues the r-expansion to 
next orders, admixtures of other irreducible representa- 
tions will appear in the formalism. Within the symmetry 
analysis, it is often argued that such terms should arise 
because the appearance of the condensate at T < Tc 
already changes the symmetry of the systemi^ The re- 
construction yields a clear quantitative estimate for such 
admixtures. It is easy to sec from Eq. ^ that the order 
parameters related to extra irreducible representations 
will be of order r"^/^ and higher, which must be neglected 
in the present analysis concerning the standard GL for- 
malism. 



IV. THREE BAND SYSTEM 

A. Eigenvectors 

As a prototype of multiband superconductors, we now 
consider a physically relevant case of a three-band sys- 
tem, the analysis of which can be done in the analytical 
form. In order to obtain the eigenvectors ^a, we write 
Eq. ([9]) as a system of linear algebraic equations 



AA 



(0) 






(0) 



0. 



(24) 



It is easy to verify that the following relations hold 



r, A(°) -9t A(") -Ti a(°) 



where 



'71 

^73 



A 723 
-4,3712 



712713, ^72 
713723, 



^2713 — 712723, 



(25) 



(26) 



and Ai is defined by Eq. ([7]). 

We now investigate the following possibilities. Let us 
first assume that 77i,f72,?/3 7^ 0. Then, from Eq. (|25p we 
immediately find that 



Ci oc 1/77^ 



(27) 



which implies that the gaps in all three bands are 
nonzero. When one of the rji^s is zero, say 771 = 0, then 
Eq. (Uni) dictates that A^°^ = A^°^ = 0, and therefore 
the condensate is formed only in one band. When two of 
the ry^'s vanish, the gap is nonzero in the corresponding 
two bands. In all these cases we deal with the nonde- 
gcnerate scenario governed by the single-component GL 



equation ((T5|) with the coefBcients given by Eq. ([161) . 
However, the eigenvector ^, whose band components ap- 
pear in Eq. (J16p . is dependent on a particular situa- 
tion. As mentioned above, for 771 , 772 , ^73 7^ we ob- 
tain Eq. ()27|) whereas for, say, rji = 772 = we have 

C^ = (1,-713/723,0). 

The case when ah 77i's are equal to zero requires a bit 
more algebra. Expressing Ai in terms of 7^ from Eq. (|26p 
and then inserting the result into Eq. ([24| . we find that 
in this case Eqs. (jM)) are reduced to a single equation 
that reads as 

7i27i3Ai°^ + 7i2723A^°^ + 7i3723Af ^ = 0. (28) 

A general solution to Eq. ((28)) can be written as 



A(0)(r)=.7?i(r)ui+792(r)H2, 



where 



Ul 




-"2 



1 

-713/723 




(29) 



(30) 



are linearly independent and 7?i,2(r) are unknown func- 
tions to be specified later. Comparing Eq. ([29]) with 
Eq. (112), we conclude that this case represents the de- 
generate scenario with A/ = 2. Equation (|29)) can be 
rewritten in terms of two orthogonal eigenvectors ^1^2 by 
applying the orthogonalization procedure to ui,2, which 
gives 



^1 = Ml, 6 = U2 - 



7?3712 



(7?2 



- 713)723 



ui. 



(31) 



One can then express the functions 'i?i_2 through -01.2 
introduced earlier as 7?i(r) = V'i(r) — 7i37i2/[(7i2 + 
7i3)723]V'2(r) and 7?2(r) = i/'2(r)- The band gaps are 
then defined by the two Landau order parameters -01 _2 
according to 



A- 



(0) 



A, 



(0) 



A 



(0) 



V'2, 
713 
712 



^1 



7?27l3 



(7?3 - 
7?37l2 



(7i3 + 712)723 



712)723 

02- 



i>2 



(32) 



Finally we note that the derivation of Eqs. (|28p - ([5^ 
assumes that 712, 713, 723 7^ 0. If some of these interband 
couplings is zero while 7?i = 772 = 773 = 0, the problem 
reduces to a trivial example of the nondegenerate case 
where some of the available bands are uncoupled. 



state with the phase frustration or the chiral solution. 
This state is of a particular interest as it breaks the time- 
reversal invariance in the system, leading to many uncon- 
ventional superconducting properties!^ Below we analyt- 
ically demonstrate the possibility of such a state in the 
three-band system within the simple variant of the model 
with strong interband couplings, i.e., ga = and gi^j > 
and iVi(O) = ^^2(0) — iV3(0). Our analytical considera- 
tion compliments numerical investigations in the recent 
Ref. [13. Such a model describes an interesting example 
of a system where the superconducting pairing is caused 
by the interband coupling and, as it is believed, may be 
relevant for pnictidesJ^ We are interested in the special 
case when different interband couplings are equal to one 
another, which may be dictated by some symmetry be- 
tween bands^ but is not necessarily limited to only this 
physical situation. Please note that many different com- 
binations of intra- and interband couplings can lead to 
a degeneracy of Tc and possible phase frustration (see, 
e.g. Ref. 122). However, in the absence of physical justifi- 
cations for such coupling matrices, we refrain from their 
analysis. 

Using the orthogonality conditions for ^q, 's and the fact 
that the band DOS's are equal, wc obtain ai2 = 021 = 
in Eq. (PO)) . Furthermore, it is obvious that the tensor 
bap-yS is symmetric with respect to the permutation of 
each pair of the indices so that it is convenient to intro- 
duce new notations 



Pi = follll, /32 = &III2 = ^1121 = ^1211 = ^2111, 

Ps ~ ^1122 = ^1212 = ^2112 = ^2121 = ^1221 = ^2211, 

/^4 = ^1222 = ^2122 = ^2212 = ^2221, /^5 = ^2222- (33) 

Then, for a homogeneous case without a magnetic field. 



Eq. (1201) yields 



ai = - PMxf - /32(27^i>2 + 7^17/^2*) 

-/33f2|7^2|^+7/'2'^)-/34|V2p^, 



'^1 



^1 



a2 



/?2|7/'lpt^-/33(2|V^lP+7A?^) 



/34(27/;i0*+ 7^1 7/^2 



k\U\ 



(34a) 



(34b) 



where we also denote a\ = an, a2 — 0,22- As usual, it 
is convenient to search for a solution to Eq. (|34l) in the 
form i^i = |7/'i| exp(z0j). Then, matching the imaginary 
parts in Eq. p4ap [or in Eq. (|34bp . which gives the same 
result], we obtain 



B. Chiral state w^ith phase frustration 

Under certain conditions the ground state of a three- 
band superconductor may develop a nontrivial phase dif- 
ference between different band gaps, referred to as the 



[/327- + 2/33 cos(<50) -f /34r-^] sin((5(/)) = 0, (35a) 

where the notations (50 = 02 — 0i and r = |7/'i|/|7/'2| are 
introduced. Matching the real parts in Eqs. p4ap and 



(I34bp yields, respectively, 



--^1 + (3/32 + Pir-^)r~^ cos((50) 
+ /33r-2 [2 + cos(25(/))] , (35b) 



IV'i 



a2 



--{I32r + 313^7--^) cos{5(j,) 
+ P3[2 + cos{2S^)] +P5r- 



(35c) 



To check the thermodynamic stability of different so- 
lutions to Eqs. pSap and ((35)) . one needs to calculate the 
free energy from the functional in Eq. ()22p . It can be 
rewritten, using the new notations, as 



F^ d 



■[\U{ 



ai+ a2T 



-'' + -MA' 



X (/3i + 4 cos{5(t))r-^ {P2 + PiT'^) 

+ 2p3r-^ [2 + cos(25(j))] + p^r-") }. (36) 

To proceed further, we substitute the chosen model 
parameters into the obtained equations. Inverting the 
coupling matrix yields ju = — l/(2g) and ji^j = l/{2g). 
Then, using Eq. ([31]) . we obtain the eigenvectors as 
^^ = (0,-1,1) and S^^ = (2,-1,-1), where ^2 is now 
multiplied by 2 for the sake of convenience of our fur- 
ther calculations. Substituting these eigenvectors into 
Eq. dUl), we find 

ai — 2a, a2 — 6a, /3i = 25, /32 = 0, 

/33 = 26, /34 = 0, /35 - 186, (37) 

where a = oi = 02 = 03 and 6 = 61 = 62 = ^3, and a^ 
and bi are given by Eq. ([S]). Finally, based on Eq. (P7|. 
we can rewrite Eq. p5ap as 



sin(2(50) = 0, 



(38) 



which yields the obvious solution for the phase difference 
S4> = 7m/2, with n being integer. 

One can identify two solution classes. The first one is 
given by (50 — 0, tt, 2tt, ... at which cos{264)) = 1. Here a 
sign difference can occur between the band components 
but there is no nontrivial phase difference. In this ease 
Eqs. psbp and p5cp are reduced to 



|^iP = |a|/[6(l + 3r-2)]. 



(39) 



The complete homogeneous solution for the band gaps is 
then given by 



aw = 




(40) 



where r ~ |'0i|/l''/'2| serves as a parameter. The corre- 
sponding free-energy density / = J^/V is obtained as 



/ = -a'/b. 



(41) 



Notice that since Eq. (|4ip does not depend on r, this 
quantity is not fixed and therefore we obtain a continuous 
family of solutions with the same free-energy density. 

The second solution class is obtained when n is odd, 
i.e., 6(j) = 7r/2, 37r/2, . . . and cos{2S(p) = —1. In this case 
Eqs. pSbp and pSep yield the system of two equations 



-5/IViP 



=6(1 



% 



-3a/|^iP=6(l + 9r-2). 



(42a) 
(42b) 



This system is solved trivially giving IV'iP = 3a/ (46) 
and r = -s/S. Then, using Eqs. P^ and taking Scj) = 
tt/2, 5tt/2, . . . and (50 = 37r/2, 7tt/2, ... we obtain two dif- 
ferent solutions as 

J2ti/3 

'(43) 

These are ehiral solutions with a nontrivial phase differ- 
ence between the band gaps. The free-energy density for 
both of them reads as 





/ = -3aV(26) 



(44) 



Comparing this with Eq. (|4ip reveals that the ehiral so- 
lution is more favorable energetically and thus the three- 
band model with strong interband couplings supports the 
formation of the ehiral state. 

This conclusion agrees with numerical simulations of 
Eq. (IH) for the three-band case,— which showed that 
the ehiral state with the phase shifts ±27r/3 is found at 
T ^ Tc only in the limit (723 — >■ 512 = 313 • The phase 
shift obtained in our work is independent of temperature, 
which differs from numerical simulations in Ref. Il8l . We 
note, however, that these simulations employed the un- 
reconstructed GL equations, where a solution does not 
account for all relevant terms of the r-expansion. A cor- 
rect temperature dependence of the phase shift must be 
calculated with the help of the extended multiband GL 
formalism that should be constructed in the spirit of the 
approach in Ref. [2^. 

As already mentioned above, the appearance of the 
ehiral state may indicate the symmetry of the model, 
reflected in the structure of the coupling matrix. In par- 
ticular, the matrix investigated in this section can be re- 
alized by choosing the bands as the pockets of the Fermi 
surface centered around X points of the Brillouin zone 
of the fee lattice, see Ref. [25|. The band gaps are then 
transformed according to a three dimensional represen- 
tation of the Oh cubic symmetry group. This representa- 
tion splits into one dimensional Ag and two-dimensional 
Eg irreducible representations. The two-dimensional rep- 
resentation Eg, that corresponds to the highest critical 
temperature, can have its basis chosen as two vectors in 
Eq. P5)) . Constructing the Landau theory from this ir- 
reducible representation, one recovers the reconstructed 



GL formalism discussed above, which additionally proves 
its validity. 

Here we stress that recasting the multicomponcnt 
GL theory ^ in terms of the basis functions of the 
symmetry-group representations does not eliminate ad- 
mixtures of different irreducible representations in the 
free-energy functionali^ However, following our analy- 
sis, such admixtures must be neglected as they exceed 
the accuracy of the GL theory, in full agreement with 
the standard Landau approach. The proper accounting 
of the admixture terms can be done only by employing 
the extended GL theory which collects all relevant higher- 
order terms in the expansion of the band gaps. 

Finally, as the chiral state is related to the degener- 
acy of a solution for Tc that can be caused by, e.g., the 
crystalline symmetry, the existence of any simple rela- 
tion between the chiral state and signs of the intcrband 
couplings ^i^j, as suggested in Refs. l20ll23 . appears to be 
very questionable at least in the GL domain. Notice that 
this conclusion is also supported by numerical investiga- 
tions of Eqs. ^ performed in Ref. U^ 



V. SUMMARY AND CONCLUSIONS 

In this work wc have derived the consistent GL theory 
from the multiband BCS Hamiltonian. The derivation 
applies a reconstruction procedure to the conventional 
Gor'kov truncation of the matrix gap equation. This 
reconstruction invokes the expansion in powers of t and 
removes incomplete contributions to band gaps of orders 
higher than r^ ", thus matching the accuracy of the gaps 
with that of the Gor'kov truncation. 

When the solution for Tr is not degenerate, we recover 
the earlier results of Refs. I3l l26l - l28l that the GL theory of a 
multiband superconductor maps onto a single-component 
GL formalism in which the spatial profiles of all band 
gaps are equivalent. However, this result is valid only 
in the standard GL domain, i.e., to the accuracy A^ oc 
T^/^. Difference between the spatial profiles of the band 
gaps appears already in the leading correction to the GL 
theory.i^Ii^a 

If the solution for Tc is degenerate, which appears due 
to a symmetry of the system, the GL theory acquires 
several order parameters. We have carried out a detailed 
analysis for the three-band system treated as a prototype 
of a multiband superconductor. For the simple three- 
band model of pnictidcs with dominant intcrband cou- 
plings, the solution for T^ is two- fold degenerate and the 
GL theory has two order parameters ipa which corre- 



spond to the two-dimensional irreducible representation 
of the relevant symmetry group, in full agreement with 
the Landau theory. We have shown that the band energy 
gaps themselves cannot be interpreted as the Landau or- 
der parameters in a multiband superconductor due to the 
Josephson-like coupling between bands. 

Our approach yields explicit expressions for the coef- 
ficients of the GL theory. Also, the formalism provides 
a solid basis for further extensions of the theory and, in 
particular, offers the correct way to account for the in- 
fiuence of other irreducible representations not inherent 
in the ordinary GL approach. 

Although it was not a purpose of our work to discuss 
the origin of the degeneracy of T^ it is worth noting 
that it does not always appear due to the crystalline 
symmetry. It can arise, e.g., in the atomically fiat su- 
perconducting nanofilms, where the size quantization of 
the perpendicular motion of electrons results in the for- 
mation of multiple single-electron subbandsi^ Such su- 
perconducting nanofilms can be regarded as effectively- 
multiband superconductors with the interaction matrijs^S 
9ij — 5(1 + ^ij/'^)/d, where d is the nanofilm thickness 
and g is the coupling constant for the material of the 
nanofilm. The structure of this matrix is similar to that 
of 5 investigated in Sec. IIV Bl and, as a result, a degener- 
ate solution for Tc also appears in this case. The devel- 
oped formalism thus provides a general link between the 
multiband BGS theory and the phenomenological Lan- 
dau model for multiband superconductors, irrespective 
of the origin of the symmetry. 

We conclude by noting that the degenerate regime 
manifests itself in several important physical conse- 
quences such as the formation of the chiral ground state 
and the appearance of different spatial length-scales of 
the band condensates, which can be observed even at 
T ^ Tc- This may result in a plethora of new phe- 
nomena, i.e., fractional vorticesjiSt fiux-carrying topolog- 
ical solitons^ and other exotic statcsi^ So far those 
phenomena have been studied using the unreconstructed 
multi-component GL model given by Eqs. ([6]), and so we 
suggest revisiting these problems in the framework of the 
true GL formalism. 
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